Meandering rivulets on a plane : a simple balance between inertia and capillarity ? 
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Experiments on streams of water flowing down a rigid substrate have been performed for various 
plate inclinations and flow rates, and we focused on the regime of stationary meanders. The outcome 
is that (i) the flow is highly hysteretic : the shape of the meanders varies with flow rate only for 
increasing flow rates, and the straight rivulet regime does not appear for decreasing flow rate, (ii) 
A simple force balance, including inertia, capillary forces, and also hysteresis of wetting, accounts 
well for the experimental instability threshold flow rate and for the final radius of curvature of the 
meanders. 

PACS numbers: 47.20.-k, 68. 08. Be, 47.60.+i, 47.20. Ma 



Meandering is ubiquitous in nature and most famil- 
iar from rivers, but can also be seen at a much smaller 
scale like on window-panes during a rainfall. Recently, 
Drenckhan et al. [l[ identified regular meandering pat- 
terns of rivulets, confined between two vertical plates, 
and charged with surfactants, mimicking undulations of 
foam Plateau borders. This observation confirms the 
generality of the pattern and suggests that a well de- 
fined mechanism is at work in rivulets, able to select a 
given meandering length scale. Understanding this mech- 
anism and the properties of the resulting meanders is by 
itself a fundamental challenge which can also have im- 
plications in other fields of physics. For instance, Bru- 
insma 2] pointed out a possible analogy with the sta- 
tistical mechanics for directed polymers inside a random 
matrix. The stability and characteristics of meandering 
rivulets may also be of importance to several industrial 
processes. In heat exchangers, changes between differ- 
ent flow regimes can cause drastic modifications of heat 
transfer [3(. Meander formation can also be an undesir- 
able feature in coating processes Q. 

For rivers @, Q , it seems asserted that erosion is the 
key mechanism, but the meandering on non-erodible sur- 
faces is still an essentially open problem. Culkin Q and 
Nakagawa [H, 0] reported significant experimental work 
on streams running down an inclined plate. Nakagawa 
identified four regimes depending on flow rate : drops, 
meanders, unstable stream (the main rivulet oscillates 
and splits into several smaller ones) and a restable stream 
(the rivulet restabilizes into a straight rivulet of variable 
width, forming a braided pattern [10|]). He stresses that 
he never obtained straight rivulets (of constant width), 
though this regime has been reported by Schmuki and 
Laso who also investigated the effects of viscosity and 
surface tension, showing that meandering is suppressed 
at high viscosities. All in all, meandering has been ob- 
served in a wide variety of configurations, but little or 
no general conclusions have been drawn. In several stud- 
ies the order parameter is the ratio of the length of the 
meanders over the length of the inclined plane (termed 



sinuosity, well kown in a geological context), a quantity 
difficult to interpret physically. 

Stability analyses of rivulets, some neglecting longitu- 
dinal flow [l2l Il3j , are mainly focused on varicose modes 
and not sinuous modes, i.e. are not dealing with the me- 
andering instability. Only two recent papers deal with 
meandering threshold 

This paper studies the shape and behavior of meander- 
ing rivulets as a function of two control parameters (flow 
rate and plate inclination), and discusses physical inter- 
pretations, focusing on the role of hysteresis. Up to now, 
there has been no quantitative study of the morphology 
of meanders on non-erodible surfaces in well controlled 
and reproducible conditions, and in addition, no com- 
parison with simple hydrodynamic models. Our paper is 
a first step in this direction. In the present paper, ex- 
periments were only performed at low viscosity (that of 
water), corresponding to common natural situations. 

Experimental set-up - Figure [T]a shows a schematic 
diagram of the experimental set-up. De-ionized water 
is injected at the top of an inclined plate (1.20 m long 
and 50 cm wide) . The substrate is a Mylar sheet (plastic 
sheets of polyethylene terephthalate (PET)), flattened on 
a rigid plate. It insures partial wetting conditions for wa- 
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FIG. 1: a) Experimental setup b) Stationary meanders for 
a = 32°. Q=1.08mL/s for the left picture and Q=1.40mL/s 
for the right one. 
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ter (advancing and receding contact angles of respectively 
9 a = 70° and 9 r = 35°) and reduces problems of static 
electricity compared to other common plastics. The tilt 
angle a of the plate can be changed at will between 0° 
and 87°. The water is collected in a tank and pumped 
back to the top of the plate by a gear pump providing 
an adjustable and constant flow rate Q, checked with a 
precision flow-meter. Pictures and movies of the experi- 
ments were taken by a digital camera placed 1 m above 
the plate and perpendicular to the latter. 

Meandering thresholds - For increasing flow rates, the 
following regimes are observed, (i) Drops. Individual 



drops periodically detach from the injector 151 ]. (ii) 
Straight rivulets. The liquid flows down forming a 
straight continuous ridge along the direction of steep- 
est descent, (iii) Meandering rivulets. Above a critical 
flow rate Q&, depending on the plate inclination, the 
straight stream is unstable. Perturbations (surface de- 
fects, injection noise, air movement, . . . ) appear as small 
bends of typical size comparable to the rivulet width w. 
These bends initially amplify laterally and downwards 
and eventually reach a stationary shape. It can take 10 
minutes to one hour for a meander to fully develop along 
the whole length of the plate, but the resulting path is 
completely stationary (Fig. [If)). Snapshots of a settled 
meander, taken every 2 minutes during 24 hours, con- 
firmed that it did not move at all throughout that time 
lapse. Using either a steady pump or a constant level 
tank, differently designed injectors and a long plate to 
check whether the behavior of the rivulet depended on 
the distance to the injector, we have verified that the me- 
andering regime existed independently of injection con- 
ditions, confirming previous findings by other groups. To 
obtain yet other injection conditions, it might be inter- 
esting to study the collision/merging of two rivulets, (iv) 
Dynamic regime. Above a second critical flow rate Q C 2, 
meanders no longer remain stable. The rivulet sweeps 
from side to side, similar to the free end of a garden 
hose [l6j |. frequently breaking up into sub-rivulets, (v) 
Restable regime. For even higher flow rates, the rivulet 
restabilizes and becomes straight again, but its width 
now varies like braids [10| . 

Fig. [2)3 displays the critical flow rates for the onset 
and disappearance of stationary meanders as a function 
of plate inclination a. The decrease of the second crit- 
ical flow rate Q C 2 with a is similar to that of Q c i, but 
we have no explanation for it yet. The dependency of 
Q c i on a can be understood from the balance of forces 
acting on the rivulet : gravity, surface tension, inertia 
and contact line pinning forces. For the lateral sta- 
bility of a straight rivulet, gravity does not intervene 
(Fig. ^M>)- Surface tension opposes the bending of the 
rivulet and, once integrated across the cross-section, can 
be seen as a line tension of the liquid rim creating 
a normal force F 7 straightening the rivulet. Taking into 
account the interfacial energies and the capillary pres- 



sure inside the rivulet, one gets F 7 = C(9)~fw/r c , where 
C{ff) is a constant (C(0) ~ 9 2 /3 in the limit of a small 
average contact angle 0), 7 the surface tension, w the 
width of the rivulet and r c the initial radius of curva- 
ture. Pinning forces are reactive (therefore stabilizing 
forces) and act normal to the contact line. Their upper 
bound is given by the advancing and receding contact 
angles: F h < F™ ax = j(cos9 r — cos6> a ). Instability will 
arise when inertia (Fi = pSv 2 jr c where p stands for the 
density of the liquid, S for the cross-section of the rivulet 
and v for the RMS velocity (~ average velocity) inside 
the rivulet) becomes stronger than both line tension and 
pinning (Fi > F 7 + _F™ ax ). The onset of meandering is 
therefore given by: 



,9k 

Sr n 



C(9)w 



+ (cos 9 r — cos 9 a ) 



(1) 



a 2 



\i 


— «-Q c1 =a 1 .(sinct)- 3 ' 5 




a,=(0.47lfl.02) mL/s R=0.96 ' 


\i 


Q 02 = a 2 


(sina)" 3 ' 5 


a2=(1 .33±0.06) mL/s R=0.93- 


V 

\ K 


J 

DYNAMIC 






MEANDERING 


B 1 1 ■ 
* I i- 




STRAIGHT 


ii- fr-i-i # 




a) «(°) 

FIG. 2: a) Meandering thresholds for increasing flow rates 
b) Forces acting on a meander, in the plane of the plate 

Without the pinning term, this balance was also sug- 
gested by Drenckhan et al. [1] for meanders in foams. 
If the pinning term can be neglected, or else if r c scales 
as w with flow rate, the first critical flow rate scales as 
pQ 2 cl / S oc jw. Let us now assume that the flow inside 
the rivulet is a Poiseuille flow (v oc w 2 g sin a / V) , where 
v denotes the kinematic viscosity), and approximate the 
cross-section of the rivulet to a disc segment with 7r/4 
contact angle (S = (tt — 2)w 2 /8). Using flow rate conser- 
vation (Q = Sv) leads to the scaling 



Qci oc (7/p) 4/5 (v/gf' a .(sin a) 



3/5 



-3/5 
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This scaling matches the experimental data very well 
(Fig [2)3). If we keep the theoretical prefactor evaluat- 
ing to about 2.3, we obtain flow rates 5 times too small. 
Rough estimates of the force terms in equ (JXJ) , using ex- 
perimental data for typical threshold conditions, show 
that pinning is not completely negligible 18]. We thus 
expect equ ([2]) to underestimate the threshold. It is satis- 
fying that the order of magnitude of the forces, and their 
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FIG. 3: Experimental data for the mean radius of curvature, wavelength and amplitude for numerous plate inclinations. 



scaling, are correctly estimated. The scaling of equ <j2j) is 
close to that predicted by Bruinsma |2j. 

If the flow rate is now decreased while in the me- 
andering regime, meanders remain remarkably stable : 
they keep their path, but become thinner, until they 
break up into drops. In particular there is no sinu- 
ous/straight rivulet transition: below Q c i the rivulets 
meander instead of becoming straight. This strong hys- 
teresis is a consequence of pinning effects. For all fea- 
tures of a meander with a radius of curvature greater 
than C(6)w/(cos6 r — cos(9 a ) ss OAw, i.e. for all but very 
small scales, contact line pinning forces (now acting in 
the opposite direction) will dominate restoring capillary 
forces (equ. {J)), even if inertial forces tend to zero. 

Shape of the stationary meanders - All the experimen- 
tal data presented here were obtained for increasing flow 
rates. For numerous meanders achieved for plate inclina- 
tions varying between 6° and 87°, we measured the mean 
radius of curvature (R c ) at the apex of the curves (mean 
taken over all curves of a given rivulet), mean wavelength 
(A), and mean amplitude (A) (see figures [1^ and [3]). All 
three parameters increase monotonically with the incli- 
nation of the plate a and the flow rate Q. 

The magnitude of inertial and capillary forces depends 
on the shape of the rivulet path. The equilibrium of grav- 
ity, inertia, capillarity and pinning forces can therefore be 
used to solve for the expected radius of curvature of sta- 
tionary meanders. Again, gravity does not matter if we 
are interested in the curvature at the vertical segments 
of the meander. The force balance at the threshold of 
depinning/motion then reads: 

(pQ 2 /S - C(0)jw) /R c = 7 (cos T - cos a ) (3) 

This is the same as equation ([1]) except that R c is now 
the final radius of curvature of the bends in the mean- 
der, which scales between 1 and 6 cm (see figure 3) and 
is much larger than the initial radius of curvature r c con- 
sidered for the threshold of the instability. An order of 
magnitude calculation [l8[ shows that the capillary con- 
tribution C (8)"fw I R c now becomes smaller than pinning, 
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FIG. 4: Comparison of the data to the model for the mean 
radius of curvature. Insert: ratio of the wavelength to the 
radius of curvature for all plate inclinations a. 

and is also smaller than inertia pQ 2 /(S.R C ). We there- 
fore neglect the capillary part and, assuming a Poiseuille 
flow, one gets the following scaling law for the radius of 
curvature: 
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This expression has been tested on our data (Fig. [4| and 
the scaling law fits the data reasonably well. The same 
scaling seems to hold for amplitude and wavelength. 

Some sets of data, for a given inclination, seem to have 
a slightly different slope from the proposed scaling law. 
A first possible explanation lies in the neglected capil- 
lary term, which would add a contribution to the radius 



of curvature in —Q ' sin - ' a, which is more important 
at low inclinations and low flow rates, and diminishes the 
final radius. A second important point is the reduction 
of the effective slope through sinuosity £>(ratio of total 
length of the rivulet to the distance between its end- 
points). The slope sin a has to be replaced by sina/iS 
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in equ ([4]), and since ^increases with Q and a @], the 
exponents of Q and a should be smaller than the ones 
in equ[U Last not least, the Poiseuille flow assumption 
might need to be corrected. Indeed, flow visualizations 
with methylene blue dye show regions with zero or even 
backward recirculating flow, already mentioned by [l9j |. 
The Reynolds number based on the rivulet size is of order 
500-1000. 

Discussion - The scenario of development of meanders 
may be summed up as follows. Stable rivulets become 
unstable when inertia dominates capillary and pinning 
forces. The destabilizing term diminishes as a perturba- 
tion grows (and so its radius), so the hysteresis of wetting 
eventually stops the growth of the perturbation. In the 
final state, pinning forces should therefore be maximally 
mobilized, with the critical advancing contact angle on 
the outside of the curve and the critical receding contact 
angle on the inside. Cross-section measurements by Nak- 
agawa and Scott [1] seem to corroborate this, but more 
quantitative measurements are planned. 

As shown in figHt, the final meanders are made up of 
circular segments. The ratio (X) / (R c ) is expected to be 
close to 4 for paths made of half-circles possibly followed 
by horizontal segments. This is confirmed by the insert 
in fig U within experimental scatter. 

In order to understand how the rivulet switches from 
one bend to the next, gravity has to be taken into ac- 
count to explain why the rivulet does not follow the cir- 
cular path back uphill. If weight is strong enough, it will 
exceed pinning at a critical orientation of the rivulet with 
respect to the slope, sin$ = 7(cos# r — cos9 a )/ pgS since 
which will define the inflection points on the path. If 
pinning always remains stronger, the rivulet can become 
horizontal without slipping : horizontal segments appear, 
in which the velocity decreases and the cross-section in- 
creases until gravity becomes dominant. Both scenari 
seem to correspond to observed patterns, but the effect 
on the amplitude of the meanders is not clear yet. 

A last point remains to be discussed : above Q c \ per- 
turbations are stabilized at large scales because the desta- 
bilizing forces decrease for smaller curvatures, but why 
do no new perturbations appear at small scales? A possi- 
ble reason might lie in the considerable drop of the mean 
flow velocity (by more than a factor 2) at the transition 
from straight-to-meandering, which has been revealed by 
injection of methylene blue dye into the rivulets. The for- 
mation of a meandering path decreases the mean slope 
seen by the rivulet, and consequently the equilibrium 
velocity. The rivulet therefore becomes stable again in 
the meandering part, and remains unstable in the parts 
which have not yet developed bends. It would be inter- 
esting to perform velocity measurements to see whether 
the second threshold Q C 2 is defined by the same critical 
local velocity as Q c \. 

Concerning the influence of viscosity on meanders, 
Schmuki and Laso [ill ] showed that viscosity damped me- 



anders. Our scaling laws are in agreement with this ob- 
servation : for higher viscosity, both the critical flow rate 
for the onset of meandering and the radius of curvature 
increase. Other experiments, varying viscosity, should be 
carried out to quantitatively determine the influence of 
viscosity on meandering. 
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